February 8, 2008 



Conservation and dynamics of blocks 10 in the cellular automaton 

rule 142 

Henryk Fuks 

Department of Mathematics, Brock University 
St. Catharines, Ontario L2S 3A1, Canada 
Toronto, Ontario M5T 3J1, Canada 
Email: hfuks@brocku.ca 



Abstract 

We investigate dynamics of the cellular automaton rule 142. This rule possesses additive invariant of 
the second order, namely it conserves the number of blocks 10. Rule 142 can be alternatively described 
as an operation on a binary string in which we simultaneously flip all symbols which have dissenting right 
neighbours. We show that the probability of having a dissenting neighbour can be computed exactly 
using the fact that the surjective rule 60 transforms rule 142 into rule 226. We also demonstrate that 
the conservation of the number of 10 blocks implies that these blocks move with speed —1 or stay in the 
same place, depending on the state of the preceding site. At the density of blocks 10 equal to 0.25, the 
rule 142 exhibits a phenomenon similar to the jamming transitions occurring in discrete models of traffic 
flow. 

1. Introduction 

Let s be a binary string of L symbols, i.e., s — soSi . . . sl-i, where s$ € {0, 1} for < i < L, L € N. We 
will say that the symbol s, has a dissenting right neighbour if Sj_i = Si =/= Si+i. By flipping a given symbol 
Si we will mean replacing it by 1 — Sj. 

Consider now the following problem: suppose that we simultaneously flip all symbols which have dissent- 
ing right neighbours, as shown in the example below. 

■•• 1 1 1 1 1 1 ••■ 

(1) 

■•• 1 1 1 1 1 1 1 ••• 

Assuming that the initial string is randomly generated, what is the probability that a given symbol has a 
dissenting right neighbour after t iterations of the aforementioned procedure? 

In order to answer this question, we will take advantage of the fact that the process described in the 
previous paragraph is actually a cellular automaton rule 142, using Wolfram's numbering scheme for ele- 
mentary cellular automata [Q . It has a property which turns out to be crucial to the solution. If one counts 
the number of pairs 10 in before and after the bit-flipping operation, it is easy to see that this number 
remains constant (similarly as the number of pairs 01). We will show that this is true for an arbitrary string, 
and we will take advantage of this fact to compute the probability of having a dissenting neighbour. 

2. Definitions 

Before we proceed, we will introduce several concepts of cellular automata theory. Let Q — {0, 1} be called 
a symbol set, and S = {0, 1} Z be called the configuration space. A block of radius r is an ordered set 
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b-rb-r+i ■ ■ - b r , where r £ N, 6, S Let r S N and let £? r denote the set of all blocks of radius r over Q. 
The number of elements of B r (denoted by card2? r ) equals 2 2r+1 . 

A mapping / : {0, l} 2r + 1 i— > {0, 1} will be called a cellular automaton rule of radius r. Alternatively, the 
function / can be considered a mapping of B r into Bq = Q = {0, 1}. 

Corresponding to / (also called a local mapping), we define a global mapping F : S — > S such that 
(F(s))i = /(s;_ r , . . . , Si, . . . , Si+ r ) for any s£ S. The composition of two rules f, g S T can be now defined 
in terms of their corresponding global mappings F and G as (F o G)(s) = F(G(s)), where s £ S. We note 
that if / £ T v and g £ T q , then / o g £ J- p + q . For example, the composition of two radius-1 mappings is a 
radius-2 mapping: 

(/ sOO-2, s-x, s , si,s 2 ) = f(g(s-2, s-x, s ),g(s-x,s 0) si), g(s , si,s 2 )). (2) 
Multiple composition will be denoted by 




n times 



A block evolution operator corresponding to / is a mapping f : B i— > B defined as follows. Let r > p > 0, 
a £ B r , f £ J- p , and let bi = f{a^p, a^ p+ i, . . . , a-i +p ) for — r + p < i < r — p. Then we define f (a) = b, 
where b £ B r - p . Note that if b £ B x then f(b) = f(b). 

In this paper, we will be concerned with trajectories of a given configuration under consecutive iterations 
of F. Denoting the initial configuration by s(0), the image of s(0) after t iterations of F will be denoted by 
s(t), i.e., 

*(t) = **(«(())), (4) 

which imples that 

s(t + l)=F(s(t)), (5) 

and hence 

Si(t + 1) = /(si_ r (t), Si_ r+ i(£), . . . , s l+r (t)). (6) 
Cellular automaton rule 142, which is the subject of this paper, has the following local function 

/(0, 0, 0) = 0, /(0, 0, 1) = 1, /(0, 1, 0) = 1, /(0, 1, 1) = 1, (7) 
/(l, 0,0) = 0, /(1,0, 1) = 0, f(l, 1,0) = 0, /(l, 1, 1) = 1, 

which can also be written in an algebraic form 

f(x Q ,xi,x 2 ) = xi + (1 - x Q )(l - xi)x 2 - x xi(l - x 2 ). (8) 



3. Conservation 

As shown in [2], rule 142 is one of the few nontrivial elementary rules which posses the second order additive 
invariant. It conserves the number of blocks 10, and this fact can be formally described as follows. Let us 
first define a function £(xo, £i) = £o(l — ^l), which takes value 1 on block 10 and value on all other blocks 
of length 2. We will call £ the density of blocks 10. Following we will say that £ is a density function of 
an additive invariant of / if 

L-l L-l 

^2 £(/( Sl ' Sl +!' ^+2),f(s t +i,s l+2 , s t+3 )) = 2J £0»i, Si+i) (9) 

i=0 i=0 

for every positive integer L and for all so, si, ■ • ■ , sl-i £ {0, 1}. In the above, and in all subsequent con- 
siderations, we are assuming that addition of all spatial indices is performed modulo L. That is, we will be 
concerned with periodic configurations, or, in other words, configurations with periodic boundary conditions 
where s(i + L) = s(i) for all integers i. 
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Figure 1: Two consecutive images of a sample configuration under rule f42, assuming periodic boundary 
conditions. Blocks 10 are labelled with letters a,b,and c. 

The right hand side of the above equation simply denotes the number of blocks 10 in the configuration 
s — (so, Si, . . . , sl_i), and the left hand side denotes the number of these blocks in the image of s under /. 
Figure ^ shows an example of a configuration s consisting of 11 sites, and its two consecutive images under 
rule 142, i.e., -£142(3) and F^ 42 (s), where F142 denotes the global function of rule 142. Periodic boundary 
conditions are assumed. The initial configuration s contains three blocks 10 labelled a, 6, and c, and one can 
clearly see that the number of blocks 10 remains constant after each application of Fi4 2 . Moreover, since 
the number of blocks 10 remains constant, we can label them with distinctive labels, which will allow us to 
keep track of individual blocks. In Figure^ such labelling can be simply obtained by enumerating blocks 10 
from left to right 1 . For example, looking again at Figure ^ we could say that block a remains in the same 
position after the first iteration, but moves to the left by one site in the second iteration. Similarly, block b 
moves by one site to the left in both iterations shown in Figure ^ 

To formalize the concept of the motion of blocks, we will first prove that rule 142 conserves the number 
of blocks 10 in an arbitrary periodic configuration. Let us note that 

s i+ i, s i+2 ), / (Sj+l, s i+2 , s i+3 )) = (10) 
/(Sj, Si+i, Si +2 )(1 - f(s i+ i,s i+2 ,s i+3 )). 

Using JHJ, the right hand side of the above equation becomes somewhat complicated, but it drastically 
simplifies if one notes that all variables Si in this equation are Boolean, and x n — x for all positive n if 
x € {0, 1}. After this simplification, one obtains 

£(/(Si, s i+2 ), f(s i+1 ,s i+2 , S i+3 )) = 

-SjSj+l + S i+ i+SiS i+1 Si + 2 - S i+ iSi + 2Si +3 . (11) 

We will now regroup terms on the right hand side 

£(/(sj, 8 i+ i, s i+2 ), /(s i+ i, s i+2 , Si+3)) = 
s t +i(l - s i+2 ) + s i+ is i+2 (l - s i+3 ) - SiS i+ i(l - s i+2 ), (12) 

and finally write the last equation as 

£(/(sj, s i+ i,s i+2 ), f(s i+ i, s i+2 ,s i+3 )) = (13) 

£( s i+l> S i+2 )) - J{s l+ i,S l+2 ,S i+3 ) + J(S{, Si+i, s z+2 ) : 

1 In general, such enumeration will not be unique, because we have to decide which block is first. However, by imposing some 
additional conditions, it is possible to obtain unique labelling. An algorithm producing such labelling has been described in |1] 
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where 

J(x ,x 1 ,x 2 ) = x x 1 (x 2 - 1). (14) 

This leads to 



L-1 L-1 



^2 £.(f(Si, s i+1 , s l+2 ), f(s t+1 , s l+2 , s l+3 )) = €(Sj+l, s i+2 ) 

L-1 L-1 
- ^2 ^( s »+li s »+2,Si +3 ) + ^ J(Si,S i+ i,Si +2 ), (15) 

8=0 1=0 

and, since Y^i=o J{ s i+ii s i+2, Si+3) = Y^,t=o J(si, Sj+i, Sj+2), the conservation condition @ follows. 

The eauation p3|l resembles continuity equation, with J playing a role of current, or flow of blocks 
10 PI- To see this, let us note that J(si(t), Si + i(t), s i+2 (t)) takes non-zero value only when Si(t),s i+ i(t), 
Si+ 2 {t) = 1, 1,0. Consider now a configuration containing block 110, surrounded by sites of undetermined 
state. Using definition of the local function for rule 142 (JJJ, we can construct partial state of the configuration 
at the next time step. Denoting by * an arbitrary value in the set {0, 1}, we have /(*, 1, 1) = 1, /(l, 1, 0) = 0, 
and /(l, 0, *) = 0, hence 

s(t) =...110... (16) 
s(t + l) =...100... 

We can clearly see that the block 10, when preceded by 1, moves by one site to the left in a single iteration. 
Similar argument could be used to demonstrate that the block 10 preceded by does not move: 

s(t) =...010... (17) 
s(t + l) =...?10..., 

where "?" denotes undetermined value. Additionally, note there is no other way to obtain the block 10 
in s(t + 1), that is, if Si(t + l)si + i(t + 1) = 10, then we must have either Si(t)si+i(t)si +2 (t) = 110 or 
Si-i(t)si(t)si+i(t) — 010. This demonstrates that indeed only blocks 110 can contribute to the current, in 
agreement with eq. (|14fl . 

4. Initial distribution 

Let us now go back to the problem stated in the introduction. In order to make the problem well posed, we 
need to define the probability distribution p from which the initial string is drawn. Since we know that the 
rule 142 conserves the number of blocks 10, it is natural to consider an initial distribution parameterized by 
the density of blocks 10. Let us define the expected value of £ at site i as 

p(i,t) = [Z( St (t),s l+1 (t))] = [ 8i (t)(l - s i+x (t))] . (18) 

Assuming that the initial distribution /1 is translation-invariant, p(i,t) will not depend on i, and we will 
therefore define p(t) = p(i,t). Furthermore, since £ is density function of a conserved quantity, p(t) is 
^-independent, so we define p = pit). 

The desired distribution parameterized by p can be obtained as follows. Let p S [0,1/2] be the target 
density of blocks 10, and let {A^}^ 1 be a collection of identical independently distributed Bernoulli random 
variables such that 

PriX, = 1) = 2p, (19) 
Pr(X 2 = 0) = 1 - 2p, (20) 

for all i <E {0, 1}. The initial configuration will be given by 

Si(0) = ( Y^Xj ] mod 2. (21) 
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Note that when Xi == 1, we obtain either a subsequence Sj(0)s,+i(Q) = 01 or Sj(0)s.j+i(0) = 10. Those 
subsequences occur at the same frequency, which accounts for the factor of 2 in eq. I|19fl • 

Let Pt(b) denote the probability of occurrence of block the b in the configuration s(t). If the density of 
blocks 10 in the initial configuration is p, then the probability of having a dissenting neighbour at time t will 
be denoted by P d i s (p, t). A site s, has a dissenting right neighbour if Si-iSiSi+i = 110 or Si_is,;Si + i = 001. 
Pdis(p,t) is therefore given by 

P di8 (p,t)=P i (110)+P t (001). (22) 

Although two block probabilities appear on the right hand side of the above definition, we will show that 
-Pdis (p, t) can be expressed in terms of a single block probability. 

As a first step, we note that the following properties are direct consequences of the definition (|2~TJl . 

Proposition 1 Let P (b) denotes the probability of occurrence of block b in the configuration drawn from 
the distribution given by $21\) . Then we have: 

(i) P (l) = P (l) = 1/2 

(ii) Po(10) = P(01) =p 

(Hi) Po(b) = Po(b), where b denotes Boolean conjugation of block b, i.e. bi = 1 — 

Rulel42 exhibits Boolean self-conjugacy, that is, replacing all zeros by ones and vice versa in the definition 
01 does not change the definition. This fact together with Proposition^iii) implies that P t (110) = Pt(001), 
hence 

P dis (p,i)=2P t (110). (23) 
Kolmogorov consistency conditions for block probabilities require that 

P t (110) + P t (lll) = P t (ll), 
P t (10)+P t (ll) = P t (l), 

hence 

P t (110) = P t (l) - P t (10) - P t (lll). (24) 
Using the fact that Pt(10) = p, we obtain 

P dis (p, t) = 2P t (l) - 2p - 2P(111). (25) 
The next result will lead to the elimination of Pt(l) from the above equation. 

Proposition 2 Let the initial configuration s(0) be drawn from the distribution given by H2J), and let s(t) 
be obtained from s(0) by iterating rule 14.2 t times, so that s(t) = P*42(s(0))- Then we have 

P t (l) = P t (0) = 1/2. (26) 

We will prove this by induction. Obviously, Po(l) = Po(0) = 1/2 by Proposition Let us assume that 
Pt(l) = 1/2 for some t. Block 1 has four preimages under fi42, and these are 001, 010, 011, 111. This leads 
to 

P i+ i(l) = P(001) + P t (010) + P(011) + P(lll). (27) 

Kolmogorov consistency conditions require that P((011)+Pt(lll) = Pt(ll), and, as remarked before, Boolean 
self-conjugacy of the rule 142 implies P t (001) = P t (110). This yields 

P t+ i(l) = P(110) + P(010) + P t (ll). (28) 

Using consistency conditions again we get 

P t+1 (l) = P t (10) + P t (ll) = P t (l), (29) 

and this, by induction hypothesis, yields Pt+i(l) = 1/2, concluding the proof. 
Proposition |21 simplifies eq. 125|) to 

P dis (p,t) = l-2p-2P t (lll). (30) 
Now the only thing left is to compute the probability of occurrence of block 111 in the configuration s(t). 
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5. Preimages 



In order to compute P t (lll), we will use some properties of preimages of the block 111. Let f 14 2(lll) be a 
set of preimages of 111 under f 142 . Then we have 

P t (Ul)= ]T P t _i(6). (31) 

Generalizing the above, we can write 

P t (lll) = Yl P °^> ( 32 ) 

6efii|(xii) 

where again f^^lll) is a set of preimages of 111 under f* 42 , i.e., under t iterations of fui. To find Pt(lll) 
using the above property, two steps are needed: first, we have to find the set of preimages of 111, and then 
to find probabilities of their occurrences in the initial distribution. Figure [2] shows three levels of preimages 
of 111. Upon inspection of this figure, two properies become apparent. 

Proposition 3 Let b be a t-step preimage of 111", that is, b £ fj^(lll). Then 

(i) The length of b is 3 + 2t; 

(ii) b ends with 111. 

The first property is an obvious consequence of the definition of f * 42 , and the second one can be easily 
proved by induction (omitted here). 

Further inspection of Figure [21 leads to the necessary and sufficient condition for a block b to be a t-step 
preimage of 111. Before stating this condition formally, we will explain it using an example. 

Consider the block b = 011100111, which is a preimage of 111 in three steps since fi42(011100111) = 
1100111, fi42 (1100111) = 00111, and f 142 (00111) = 111. Let us now assume that we start with a "capital" of 
1. We will move along the string b = bgbi . . . bg starting from i = 6 and moving in the direction of decreasing 
i. Every time we see that is different from bi, we decrease out "capital" by 1. If bi-i = bi, we increase 
our "capital" by 1. We stop at i = 1. 

Clearly, it is possible to traverse b = 011100111 following this procedure without making the capital 
negative. It turns out that this is a general property of preimages of 111. If b is a preimage of 111, then it is 
possible to traverse it keeping the capital non-negative. If b is not a preimage of 111, the capital will become 
negative at some point. A more formal statement of this property is as follows. 

Proposition 4 Let t be a non-negative integer, and let b = b bi . . . &2t+2 be a binary string of length 3 + 2t 
ending with 111. Define x to be a function of two variables such that x(% v) = 1 if u — v, and x(^, v) = — 1 
otherwise. The string b is a preimage of 111 under f* 42 if and only if the inequality 

k 

l + X)x(&3t-i-i,&3t-i)>0 (33) 

i=0 

is satisfied for all k = 0, 1, . . . , It — 1. 

Instead of proving this proposition directly, we will show that it can be derived from a similar result 
previously obtained for a related cellular automaton rule. 

6. Rule 226 

In |S] it has been observed that 

/60 /l42 = /226 ° /60, (34) 



G 




Figure 2: Tree of preimages of the block 111 under rule 142. 
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Figure 3: Spatiotemporal patterns generated by the rule 142 and their transformations under the rule 60. 
Densities of blocks 10 are (a) p — 0.23, (b) p = 0.25, and (c) p — 0.26. All patterns show 100 consecutive 
iterations of a lattice of 100 sites with periodic boundary conditions. Black squares represent 1, white spaces 
represent 0. 



This means that there exists a local mapping (rule 60) which transforms rule 142 into rule 226. The above 
correspondence between rules 142 and 226 has been illustrated in Figure |3J which shows spatiotemporal 
patters of rule 142 (left panel) as well as images of these patterns under the rule 60 (right panel). The 
patterns in the right panel are in fact identical to spatiotemporal patterns which one would obtain by 
iterating rule 226, providing that the initial configuration in the right panel has been obtained by applying 
rule 60 to the corresponding initial configuration from the left panel. 

Rule 226 and its image under spatial reflection, rule 184, are the only non-trivial elementary number- 
conserving rules, and many results regarding their dynamics have been established [B1I71 IH1 1^1 fTUl fTTl fT^l fTl fO | . 
For our purpose, one such result will be particularly useful. 

Proposition 5 Under the rule 226, t step preimages of 00 have the following properties: 
(i) In each preimage, the number of zeros exceeds the number of ones. 



where 



feo{xo,Xi,x-2.)=xo + Xi mod 2 

f226(xo,Xl,X 2 ) = XqXy - X±X 2 + X 2 - 



(35) 
(36) 
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(ii) The block a^a\ 



o-2t+i is an t-step preimage of 00 if and only if it ends with two zeros and 



1 + (*»+!-<) >« (37) 

for every 2 < k < 2t + 1, where £(0) = 1, £(1) = -1. 
(Hi) The number of t-step preimages of 00 containing exactly uq zeros and n\ ones is equal to 



no + m \ Hi 
where hq + n\ = It + 2. 

Proof of this result can be found in ^3], with further generalization in [TT]. The proof is based on the fact 
that the enumeration of preimages of 00 under rule 226 (or 184) is equivalent to the problem of enumeration 
of planar lattice paths between two points, subject to some constraining conditions. This path enumeration 
problem can then be solved using combinatorial methods. 
We will now relate preimages of rules 226 and 142. 

Proposition 6 The number oft step preimages of 111 under rule 14-2 is equal to the number of t step 
preimages of 00 under rule 226. 

We will explicitly construct a bijection T between f{^|(lll) and ^§(00). Let a; be a block of length m, 
60&1, ■ ■ ■ j frm-i- We define 

[T(a;)]j = Xi + x i+ i mod 2 (39) 

for i = 0, 1, . . . , m — 2. T(x) is therefore a block of length m — 1. Since T is a block evolution operator of 
rule 60, relationship similar to (|34J) must hold, i.e. 

T o f 142 = f 226 ° T. (40) 

Let us now assume that b is a t-step preimage of 111 under fi42. This means that ffxzify — HI- Since 
T(lll) = 00, we have T(f 1 ' 42 (6)) = 00. Using gDJ we obtain f| 26 (T(&)) = 00. This means that if b is a t-step 
preimage of 111 under fi4 2 , then T(b) is a t step preimage of 00 underlOO rule 226. 

Now let us consider a transformation inverse to T . In general, T is not invertible, but if restricted to the 
set of preimages of 111 under fi4 2 , it becomes invertible. 

For an arbitrary block y, there exist two different blocks x such that T(x) = y, and one can show that 
these two blocks are related by Boolean conjugacy. For example, we have T(lll) = 00 and T(000) = 00. 
We have to define T _1 such that this ambiguity is removed. This can be done as follows. Let a be a t-step 
preimage of 00 under rule 226, a = aoai . . . o^t+i- We define 



[T-^W = 1 

[T _1 (a)] 2t+ i = 1 + a 2(+ i mod 2 

[T _1 (a)] 2t = 1 + a 2t + a 2t+ i mod 2 

[T _1 (a)]o = 1 + a a + ai + . . . + a 2 t+i mod 2, 



or in a general form 



2t+l 

,l+y^a 7 - mod 2 if i = 0, 1, . . . , 2t + 1, 
\T '(,)], { ^ 3 ■ (41) 

1 if i = 2t + 2 

One can easily show that the above transformation is indeed an inverse of T, and in addition we guarantee 
that when a ends with two zeros, T _1 (a) ends with three ones, are required for an t-step preimage of 111 
under rule 142. 
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Now, if a is a i-step preimage of 00 under rule 226, we have f^g (a) = 00, hence f226C^(^ -1 ( a ))) = 00, 
and by eq. gQJ we obtain T{i[ i2 {T- 1 {a))) = 00. The last equation implies that fl 42 (T' 1 (a)) = 111, which 
means that T _1 (a) is an i-step preimage of 111 under rule 142, as required. □ 

Proposition Q] follows form the above result and Proposition EJii). 

7. Probability of occurrence of 111 

The bijective transformation T constructed in in the proof of Proposition has a property which will be 
useful in computing P t (lll). Let us call the block x = xqX\ a matching pair if xq — x\, and a mismatched 
pair if xq ^ x\. If a = T(b), then the number of matching pairs in b is equal to the number of zeros in a, 
while the number of mismatched pairs in b is equal to the number of ones in a. This fact, together with 
Proposition immediately leads to the conclusion that under the rule 142, the number of i-step preimages 
of 111 with exactly no matching pairs and n% mismatched pairs is equal to 

no — n% /no + n\ 
no + ni \ ?ii 

where no + n\ = 2t + 2. Probability of occurrence of a matching pair in a in the initial configuration drawn 
from the initial distribution is 2p, and the mismatched pair is 1 — 2p. Therefore, the probability of occurrence 
of a block with prescribed sequence of matching and mismatched pairs such that it has exactly no matching 
pairs and n\ mismatched pairs is equal to (2p) ni (l — 2p) n °. This implies that the probability that a block 
of length 2t + 3, randomly selected from the distribution (|21J) . is a t-step preimage of 111 with exactly no 
matching pairs is equal to 

^(^y r(1 - 2 ,r. (43 , 

The factor 1/2 in front comes from the fact that there are always two strings with a given sequence of pairs 
(related by Boolean conjugacy), but only one of them is a preimage of 111. 

The smallest possible number of matching pairs in a £-stcp preimage of 111 is t + 2 (recall that the number 
of matching pairs must exceed the number of mismatched pairs), while the maximum possible number is 
2t + 2 (all zeros). Summing (|43|) over no we obtain 

Pt(m) _ - .■.-(» + »-.) ( * + _« J (2p) „«-„„ (1 _ 2p) , m 

Introducing a new summation index j = no — (i + 1) we get 

laiD-gs^^^M^-'a-vr^, (45) 

and as a result, the probability i|30fl becomes 

P dis (p, *) = 1-2P-E(ii(t + t-j) P')* 1 "^ 1 ~ 2 P) t+1+J > ( 46 ) 

where p G [0, 1/2]. 

8. Equilibrium probability 

We will now show how to obtain the equilibrium probability, i.e., lim^oo Pdis(p> t). In order to find the limit 
limt_>oo -Pt (HI) we can write eq. B45JI in the form 

t+i 

p t (ui) = 2TT2 b{t + 1 - J ' 2{t + 1); 2/9) ' (47) 

3=1 
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where 

b(k,n,p)=(^jp k (l- P r- k (48) 

is the distribution function of the binomial distribution. Using de Moivre-Laplace limit theorem, binomial 
distribution for large n can be approximated by the normal distribution 

b(k, n, p) ~ 1 exp ~ (fc .; np) . 2 . (49) 

^/2itnp(\-p) 2np(l-p) 

To simplify notation, let us define T = t + 1. Now, using (|49Jl to approximate 6(T — j, 2T, 2p) in (|47(1 . and 
approximating sum by an integral, we obtain 

, x f T x 1 -(T-a;-4Tp) 2 , , 

Pt(Ul)= / — == =exp ^— - - " dx. (50) 



i 2T ^8vrZ>(l - 2p) 8Tp(l - 2p) 



Integration yields 

#(111) = 



p(l-2p)f /-(l -T + 4pT) 2 \ / -2pT 



2^T 1 6XP ^ 87X1-2/,) J CXP V2(l-2p) 
-fl - 4 ) /erf f 4pT \ _ ^ { ( 1 — T + 4pT 



where erf (x) denotes the error function 



2 f x 
erf(x) = — / e~ t2 dt. 
Jo 



(51) 



The first term in the above equation (involving two exponentials) tends to with T — > oo. Moreover, since 
lim-r^oc erf (a;) = 1, we obtain 



lim P 4 (lll) - ~(1 - 4p) J 1 - lim erf f 1 T + 4pT - ] 



Now, noting that 

/ 1 - T + ApT \ _ f 1, if4p>l, 



lim erf , = <^ \ p ' (52) 

t^oo I i/8p(l - 2/5)7 I I -1, otherwise, 



we obtain 

1/2 -2p if p< 1/4, 



lim P t (lll) = ^ I/ p T • (53) 
t^oo otherwise. 



The final expression for the equilibrium probability becomes 



9. Current 

The equilibrium probability calculated in the the previous section exhibits a singularity at p — 1/4. This 
singularity is of a similar nature as the jamming transition observed in CA rules 184, 226, and related models. 

Recall that in section 3 we defined the current J (eq. I14f> . The expected value of the current is i- 
independent, so we can define the expected current as 

j(p,t) = E„[J( Si {t),s i+1 (t),s i+2 (t))] = [ Si {t)s i+1 (t)(s i+2 {t) - 1)] . (55) 
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The graph of j(p, oo) as a function of p is known as fundamental diagram. Using the notion of block 
probabilities we can rewrite l|55|) in an alternative form as 

j{p,t) = -Pt(llO), (56) 

and using l|23[l 

i(p,*) = -^dfa(p,t). (57) 

The probability of having a dissenting neighbour, as we can see, is proportional to the expected current. 
Since the current J represents the flow of blocks 10, the expected current must be equal to 

j{p,t) = pv{p,t), (58) 

where v(p,t) is the expected velocity of a block 10 at time t. Using l|54T this velocity is given by 

We can see that for densities of blocks 10 smaller than 1/4, the average velocity remains constant and equal 
to —1, which means that all blocks are moving to the left. At p = 1/4 a jamming transition occurs, and when 
p increases beyond 1/4, more and more blocks are stopped. This phenomenon is very similar to jamming 
transitions in discrete models of traffic flow, which have been extensively studied in recent years (^3 and 
references therein). 

10. Conclusions 

We investigated dynamics of the cellular automaton rule 142. It can be transformed into rule 226 by a 
surjective transformation, which turns out to be invertible if restricted to preimages of 111. This transfor- 
mation allows to compute the probability of having a dissenting neighbour, which, in turn, allows to compute 
the expected current of blocks 10. Rule 142 exhibits jamming transition similar to transitions occurring in 
discrete models of traffic flow. 

It is worth mentioning that there are other CA rules conserving the number of blocks 10 which also 
exhibit singularities of fundamental diagrams, for example rules 35 and 14, as reported in [2]. For these 
rules, however, there exist no transformation relating them to other rules with singularities, thus the method 
presented in this paper cannot be easily applied. Nevertheless, the nature of singularities in these rules 
appears to be the same, thus some relationship between them and rules 184/226 may exist. This problem is 
currently under investigation and will be reported elsewhere. 
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